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Model quantum Hall states including Laughlin, Moore- Read and Read-Rezayi states are general- 
ized into appropriate anisotropic form. The generalized states are exact zero-energy eigenstates of 
corresponding anisotropic two- or multi-body Hamiltonians, and explicitly illustrate the existence 
of geometric degrees of in the fractional quantum Hall effect. These generalized model quantum 
Hall states can provide a good description of the quantum Hall system with anisotropic interactions. 
Some numeric results of these anisotropic quantum Hall states are also presented. 
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I. INTRODUCTION 

Variational wave functions play a fundamental role in 
our understanding of fractional quantum Hall (FQH) ef- 
fect, with the Laughlin wave function^ being the most 
prominent example. New classes of quantum Hall trial 
wave functions have been proposed as correlation func- 
tions in various conformal field theories by general- 
izing the clustering properties of the wave functions 
or as Jack polynomials with a negative parameter and a 
matching root configuration.^ For a long time it has been 
understood that the Laughlin wave function, as well as 
other FQH trial wave functions, contains no variational 
parameter. This understanding is the consequence of 
our interest in searching for topologically distinct quan- 
tum Hall wave functions; geometry, as oppose to topol- 
ogy, is believed to be redundant. Recently one of US'" 
pointed out that such topological description of FQH 
wave functions is not complete; in his geometrical de- 
scription the original Laughlin wave function is simply 
a member of a family of Laughlin states, parameterized 
by a hidden (continuous) geometrical degree of freedom. 
The family of the Laughlin states, with the geometrical 
factor as a variational parameter, should provide a bet- 
ter description of the FQH effect in the presence of ei- 
ther anisotropic effective mass or anisotropic interaction, 
which are present in real materials. 

The anisotropic FQH states are of present interest, 
both theoretically^ii^ and experimentally. In Ref. [ol, 
some of us studied the quantum Hall effects in a fast ro- 
tating quasi-two-dimensional gas of polarized fermionic 
dipoles. The fast rotation is equivalent to the high mag- 
netic field according to the Larmor theorem. And since 
p-wave interaction for the polarized fermions is typically 
very small unless in the resonance regime, the only sig- 
nificant interaction is the dipole-dipole interaction that 
could be tuned by adjusting the applied electric and mag- 
netic fields. By tuning the direction of the dipole moment 
with respect to the z axis, the dipole-dipole interaction 
becomes anisotropic in the x-y plane. Specifically, the 



system contain two parameters: the confining strength 
determined by the rotation frequency and the tilt angle 
determined by the applied field. Thus such systems have 
highly tunable anisotropic interactions, and are ideal for 
studies of anisotropic FQH states. Results in Ref. [ol 
clearly indicate the inadequacy of the known variational 
wave functions in the description of such states. 

As a matter of fact, consideration of possible 
anisotropic FQH states has a long history^ — with vari- 
ational wave functions constructed that are straightfor- 
ward modifications of the original Laughlin wave func- 
tion J^ii^ Unlike the original Laughlin wave function and 
the states to be discussed here however, these earlier 
anisotropic FQH states are not exact eigenstates of spe- 
cial 2-body Hamiltonians. Interest in such states was par- 
tially motivated by the observation of compressible states 
with strong anisotropic transport properties in high Lan- 
dau levels (LLs)J^ii^ In all these states, the rotation 
symmetry is broken spontaneously. Very recently, an 
anisotropic FQH state has been observed at u = 7/3 
In this case the rotation symmetry is broken explicitly by 
an in-plane magnetic field, whose direction dictates the 
transport anisotropy. Generalizations of the FQH states 
to lattice models with zero net magnetic field and full 
lattice translation symmetry is also called for to describe 
fractional quantum anomalous Hall states and fractional 
topological insulatorsJ^ 

Ref. pointed out the existence of a family of Laugh- 
lin states that are zero energy ground states of a fam- 
ily of Hamiltonians consisting of projection operators, 
which depend on a parameter called guiding center met- 
ric. However the states themselves were not constructed 
explicitly, and only some of their qualitative properties 
were mentioned briefly. The main motivation of this work 
is to explore the construction of a family of wave func- 
tions in closed form with numerical comparisons to facil- 
itate the study of the geometrical aspects of anisotropic 
FQH states as a result of anisotropic interaction in the 
planer geometry. 

The rest of the paper is organized as follows. In Sec.HIl 
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we introduce the anisotropic LL basis states using the 
unimodular transformation. In Sec. IIIII we focus on the 
recipe to generate the anisotropic many-particle states 
due to the unimodular transformation. Section. [TVl covers 
the numerical studies on various properties of the family 
of Laughlin states, including their density profile, pair 
correlation function, and variational energy. We summa- 
rize the paper in Sec. |Vl 



II. ANISOTROPIC LANDAU LEVEL BASIS 
STATES 

The original Laughlin wave function was most easily 
written down on a disc, using the symmetric gauge in 
which the single particle basis states are angular momen- 
tum eigenstates. The key to explicitly constructing the 
anisotropic Laughlin states proposed in Ref. 7 (on a disc) 
is the usage of a set of anisotropic LL basis states. We 
will, however, continue to use the symmetric gauge, in 
which the lowest LL (LLL) wave functions are holomor- 
phic. The same set of basis states appeared earlier in 
the consideration of deformation of shape of quantum 
Hall liquids in the context of Hall viscosityi^, although 
the wave functions were not explicitly given. In the fol- 
lowing we generate them explicitly. In addition, we also 
use them to construct the corresponding integer quantum 
Hall wave function, as a warm-up for their later applica- 
tion to FQH states. 

Let us start from the eigenvalues problem of a two- 
dimensional charged particle subjected into a uniform 
magnetic field B = Bz. The one-body Hamiltonian is 
given by 



(1) 



where we used the Einstein convention, ruah is the cy- 
clotron effective mass tensor, e > is the charge of the 
particle, A is the vector potential of the uniform mag- 
netic field and the dynamical momentum tt satisfies 



[7ra,7r5] = ieahh ^ 



(2) 



with eah = e the 2D antisymmetric Levi-Civita sym- 
bol and I = ^Jh/{eE) the magnetic length. In terms 
of a complex vector u?, the effective mass tensor can be 
written as 

rriah = m {uluh^ ulua), (3) 
{m-^f^ = m-^KV+a;^*a;^). (4) 

we note that the assumption that it is isotropic in the 
(x, y) Cartesian coordinate system means 



(5) 



The one-body Hamiltonian can now be expressed as 

Ho = lhujc{b'^b^bb'^) . (6) 



where hUc = /ml'^ is the cyclotron energy and the Lan- 
dau orbit ladder operators are given by 

b = h-^£ (u; • tt) , 6^ = h-^£ (cj* • tt) , (7) 

with [b^b^] = 1. By introducing the "guiding center" 
coordinates R, 



which satisfy 



ab n2 



[R\TTb]=0, 



(8) 



(9) 



we could define a second set of "guiding center" ladder 
operators, a and that commute with b and b^ (and 
hence with Hq) 

a = (cj* . i^) , = r ^ (cc? • i^) , [a, a^] = 1. (10) 

To proceed further, let us choose the vector potential 
of the uniform magnetic field in the symmetric gauge. 



A=^Bxr = ^B{-y,x), 
and define the complex coordinate. 



(11) 



(12) 



Let us note that for the case of isotropic mass the complex 
coordinate is explicitly given by 



X -\- iy 



(13) 



which is not the most standard definition in the liter- 
ature. Then we could express the Landau orbit ladder 
operators as (neglecting the trivial —i and i). 



b=\z^d. 



-z" - 3 



(14) 



where dzf{z,z*) is the partial derivative df /dz\z*^ etc, 
and the "guiding center" ladder operators as 



a = ^z* + dz, = \z — dz* . 



(15) 



The Hamiltonian Hq has a rotational symmetry gen- 
erated by 



Lo = a^a-b^b, [Lo,Ho] = 0, 



(16) 



where 



[Lo, x] = iy, [Lo, y] = -ix, (17) 

[Lq.Px] = ipy, [Lo.Py] = -ipx' (18) 

Simultaneous diagonalization of Hq and Lq allows a 
complete orthonormal basis of one particle states |V^nm) 
which used to be constructed as 

Holi^nm) = (n + ^)hjc\i^n7n), (19) 

^olV^nm) = (m - n)|V^nm), (20) 

\i;nm) = ^ 'J-,' IVOO), (21) 

Vm'.ni 

alV'oo) = ^>IV'oo) = 0. (22) 
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This is, of course, the commonly used basis states asso- 
ciated with the symmetric gauge. 

However, it is possible to construct a set of closely 
related but different set of basis states. First let us write 



Lo = L + Z, (23) 
L = ^(6^6 + 66^), Ho = hujcL, (24) 
L = Ka'^a^aa'^). (25) 

The conventional basis is then just the set of the simul- 
taneous eigenstates of Hq and Z, 



m^nm) = (m + ^)|^nm)- 

We may write L = L{g^), where 

L{9) = ^gabR''R\ 



(26) 



(27) 



and Qah is a positive-definite Euclidean metric tensor with 
det^ = 1. Note that L{g^) uses the "Galilean metric" 
derived from the effective mass tensor. 



(28) 



However, as has been stressed by one of us^i^, that be- 
fore specifying the two-body interactions, there is no fun- 
damental reason to choose a guiding-center basis that is 
adapted to the shape of the Landau orbits. There is a 
more general family of bases, parameterized by a complex 
number 7 with I7I < 1, 

u;a = {l-\l?)-"\u:a+l*<l (29) 
that determines a unimodular metric gah{l)^ 

1 /(l + 7)(l+7*) ^(7-7*) 
1-7*7V ^(7-7*) (l-7)(l-7*) 



according to 



(30) 



and then gah{l) = CjlCdb-\-Cja^b' then construct a 

basis of eigenstates \ipnm{'y)) of Hq and ^(7) = £(^(7)), 
which is accomplished by defining a new set of "guiding 
center" ladder operators, also considered in Ref. 



R 



and explicitly as 



V2e ' 



w R 



1 7 



(31) 



(32) 



This is a Bogoliubov (or "squeezing") transformation 
that preserves the commutation relation [a^,ajj,] = 1. 
Then the new simultaneous eigenstates of Hq and £0(7) 
are given by 



|V^nm(7)) = 



#00(7)), (33) 



Vnlml 



«7lV^oo(7)) = ^IV^oo(7)) = 0- 



(34) 



Note that \ipnm{^)) is the conventional basis \ipnm) in 
Eq. (ini). 

Since the model quantum Hall states are often defined 
in the lowest Landau level (LLL) (n = 0), it is useful 
to examine the wave functions V^om(7,''") = (^|'0Om(7))- 
First of all, from the LLL property bi/jomil ^ 'f^) = 0, we 
have 



i^omil^r) = /(2;)V^oo(0,r) 



(35) 



with 



^oo(0,r) = (27r)-i/2e-l^l'/2. (36) 
Second from a^ipoo{'y^r) = 0, we immediately obtain 

V'oo(7,r-) = Ai/2e-57^^^oo(o,r-) (37) 

with 

(38) 



Next let us note that tljomij^r)^ 

[^z-5,*+7*(^^* + S.)]^ 



could be simplified by virtue of the property that 
6V^oo(0,r) = 0, or 

5.* [V^oo(0,r)/(z)] = (-|z)V^oo(0,r)/(z) (39) 

and aipoo{0^r) = 0, or 

d. [i^ooiO, r)f{z)] = Voo(0, r)(9, - ^z*)f{z) (40) 

and 



Thus we obtain 



'm 



■{d.-^z)f{z). 



(41) 



^-0^(7,^) = i^m{l,r)^{z + z'^d,ri (42) 
^o^{^,r)^W^{z,zl) (43) 



with 



!*/>? 



(44) 



and Wmiz^Zo) are the polynomials defined by Wq{z^Zq) 
= 1, and the recursion relations, 

(45) 

W^{z,zl) = d,Wm+i{z,zl)/{m+l). (46) 

The first few of Wm{z, Zq) are given by 

Woiz,z^) = 1, 

Wiiz,z^) = z, 

W2{z,zi) = z^ + zl 

Ws{z,zl) = z^ + Zzlz, 
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and the general solutions are given by 

_ 7.!/o^2\/c r (~^) 



k\{2zlYLl-^\-z'/2zt,), (47) 
k\{2zlfzL^}\-z^/2zl), (48) 



where L^^\x) are generalized Laguerre polynomials de- 
fined by 



r(a + /c + l) 



{-xY 



Via ^ i ^ l)T{k - i ^ I) i\ 



(49) 



These expression do satisfy the reclusion relation, which 
could be proved by using the following useful identities 
for 4")(x), 



(50) 
(51) 



xL^i'^x) = {k^a)L^^},{x)-kLf\x). (52) 

Note that for \z\^ > \zll Wm(^,^o) ^ 

The densities |V^om(77 ''")P defined by the single-particle 
orbit als have some noteworthy features. First, the den- 
sity profile is anisotropic for 7 7^ 0. Take a simple exam- 
ple like V^oo(|7l,^), 



|V^oo(|7l,^)P ^ exp 



(l + |7|)x2 + (l-|7l)^' 



2P 



For this Gaussian wave packet, the width along ?/-axis is 
larger than the width along x-axis, so non-zero I7I causes 
stretching the wave function along some direction. Sec- 
ond, |'^om(77 '^)\^ has m zeros as can be seen from Fig. [T] 
in which we plot the density profiles |'?/;om(l/2, r)p for 
m = 0,1,2,3. These zeros are roots of Wm{z,ZQ), and 
they are aligned along the stretched direction. In the 
limit of 7 ^ 0, they collapse to a multiple root at the 
origin. Third, the stretching direction is determined by 
the phase of the complex number 7, arg(7), which is 
easy to understand by noting that Ws{z^Zq) = gives 
rise to three roots, and izV^zo- In Fig. [21 we plot 
|V^03(e^^/2,r)p with (j) = 0,7r/2,7r and 37r/2. Thus we 
realize that the density profile \il^om{l^'f^)\'^ is equivalent 
to rotating clockwise the density profile |'?/;om(|7|5 ''")P by 
arg(7)/2. 

To proceed further, let us define here the non- 
interacting basis states to prepare for the many-body 
problem. For bosonic system, the non-interacting basis 
state in the generalized LLL is given by 



A^^(7) (X perm [V^om(7,^)] , 



(53) 



where /i = [/ii, /i2, • • • , Mat] is a non-decreasing parti- 
tion and perm[?/^o/i(7, ''")] is the permanent of the square 
matrix whose matrix elements are '0/ii(7,Tj) (i,j = 
1, . . . , N). For the fermionic system, the non-interacting 
basis state is defined by a Slater determinant. 



=det [ipofiil^r)]. 



(54) 




FIG. 1: (Golor online) The density profiles of the generalized 
LLL wave functions tl^om{j — 1/2, r) with m = (a), m—1 
(b), m = 2 (c) and m — 3 (d). 





FIG. 2: (Color online) The density profiles of '003(7,''") with 
7 = e^'^/2 for = (a), = 7r/2 (b), = tt (c) and = 37r/2 
(d). Note that the zeros are localized at 2; = and zbiVS^o- 



where /i is a decreasing partition and det denotes the 
matrix determinant. 

As a special example, the Slater determinant ^[^0(7) 
with the partition /io = [A^ — 1, A/" — 2, . . . , 1, 0] is the 
generalized A'-particle IQH state, whose (unnormalized) 
expression could be simplified as 



(55) 



as lower-degree polynomials must vanish due to the Pauli 
exclusion principle (this is generally not the case for 
other partitions, which is the key to understanding the 
anisotropic FQH states). 

The generalized IQH state can also be viewed as a co- 
herent superposition of the isotropic IQH state (or sim- 
ple Vandermonde determinant state) and its edge states 
with angular momentum differences being multiples of 
2h. This is achieved by expanding the exponential factor 
in Eq. ([55]) . From the viewpoint of the superposition, 
this is consistent with the numerical result in Ref. 0. 



III. ANISOTROPIC MANY-PARTICLE STATES 

In this section, we first consider the two-particle prob- 
lem and then simplify the interaction Hamiltonian pro- 
jected into the LLL in Subsec. IIII A[ More importantly, 
we present the recipe to find the anisotropic counter- 
part of a LLL state and elaborate it in Subsec. IIIIBi 
We specifically consider the anisotropic counterpart of 
the prominent FQH states such as the Laughlin state in 
Subsec. [mCl 
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A. Two-particle problem and the projected 
two-body interaction Hamiltonian 

If the interaction is dominated by the energy gap huOc 
between the lowest and second Landau level, the inter- 
action Hamiltonian could be projected into the lowest 
Landau level. In terms of non-commuting "guiding cen- 
ter" coordinates Ri, the projected two-body interaction 
Hamiltonian we consider here is given by^ 



J^int({Kx},7) = ^Kx/lm(7), 



(56) 



where Vm is the anisotropic Haldane's pseudopotential 
and the projection operator hm{j) is given by 

i<j 

vUq,g) = LU9''\aqd^)eM-\9"\aqb(!'), (58) 

where Lm{x) = Lm~^\x) are the Laguerre polynomials. 
In addition, the total angular momentum associated with 
the guiding center is given by 



(59) 



For two-particle, the LLL states that simultaneously 
diagonalize £(7) and the interaction Hint{{Vk}^j) are 
completely fixed by symmetry. To clarify this, let us 
note that 



ai27«127 + ^127^127 + 1 



(60) 



with ai27 = (ai^ — a2j)/ and = (ai^ + cl2j)/ V^, 
and the only term containing operator in the interaction 
Hamiltonian, Ri — R2 , are given by 



Ui — ^2 



(1 -7*)«127 + (1 -7)«l27 



^1-^2 = [(1+7)^27 -(1+7*K27 



and contains only ai27 and a\2^. Therefore, we could 
expand £(7) and i^int({Vfe}, 7) in the orthonormal basis 



\n/r \ (^127)^(<^127)'^ |^^\ 

vM!m! 
ai^|00)^ = a27|00)^ = 0, 



and find that 



H,^t{{Vk},-f)\Mm)^ = Vm\Mm)^, 
L{-f)\Mm)^ = (M + m + l)\Mm)^. 



(61) 
(62) 



(63) 
(64) 



Thus the projected operator hmij) for many-particle 
system could be rewritten as 

= E ( E \Mm)^^Mm\\ {i,j) (65) 

i<j \ M ) 

and in second quantized form as 



M 1711171277131714 



(66) 



with |mi,m2) = |V^o,mi(7)) ^ IV^o,m2(7)) and gl^ being 
the creation operator that creates a particle in state 

|V^0m(7))- 

For the rotationally-invariant case, 7 = 0, 



^Mm(0,ri,r2) 



(^1 + ^2)^^ {zi-z^y 



n^oo(0,r,) 



and the translation into the Heisenberg state is 

|M™)7 = H±ii:^(^i^|00)7. (67) 

Therefore we find the anisotropic deformation is now eas- 
ily achieved by the simple replacement 



4 4^, 



(68) 



and the vacuum state |^o(0)) |^o(7)) that satisfy 
aijl'^oil)) =0. In the next subsection, we shall uti- 
lize this idea to consider the anisotropic deformation of 
a general LLL wave function. 



B. Anisotropic many-particle states 

The general form of a rotationally-invariant A/'-particle 
LLL wave function in the symmetric gauge is 

^(n, . . . , riv) = F{zi, ...,zn)Y[ ^oo{ri), (69) 



where F(zi, . . . , zat) is a homogeneous multivariate poly- 
nomial with F{Xzi, . . . Xzn) = F{zi^ . . . , zn) and 
1(0) I ^) = (M H- ^A^)|^). The translation to the Heisen- 
berg picture yields 

|*(7 = 0)) = F(al,...,ajv)|*o(7 = 0)), 
ai|*o(7 = 0)) = 6i|*o(7 = 0)) =0. (70) 

The deformation of this Heisenberg state becomes 

|*(7)) = F(a;,i,..., 4^^)1*0(7)), 
07,^1*0(7)) = &.|*o(7))=0. (71) 
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The Schrodinger wave function for this state |^(7)) wih 
have the form, 

^(7, ri,...,rN) = F^{zi,...zn)Y[ ^ooin). (72) 

i 

We now need the formal construction of F^{{zi})^ 
which is stih holomorphic, but no longer a polynomial. 
Using the homogeneity of the polynomial F, the wave 
function is given by 

A-^F({iz, - a,. + 7*(^< + d,,m n Xh-i-'^hooiri). 

i 

Using the property aiipoo{ri) = 0, or 

dMr)f{z, zn = ^oo{r){d, - ^^*)/(^, ^*), (73) 
and bill; 00 (vi) = 0, or 

dz*i^oo{r)f{z, z*) = iJoo{r){dz^ - \z)f{z, z*), (74) 

the holomorphic part F^{{zi}) of the wave function be- 
comes 

i 

Using the following equality, 

d. (e-5^^?/(^)) = e-^-^"- {d, - 7^) f{z), (75) 
turns into, 

which is easy to understand by noting Eq. (|42]) . Using 
the homogeneity of F, we obtain the final result 

F^({A-i^,})= (n,A^e-*^^.')F({^, + ^o'a,J)l. (76) 

where A and Zq are fixed by 7 through Eqs. (|38]) and 
(|44|) . This expression no longer requires homogeneity of 
F{{zi})^ and is quite general for the deformation of the 
holomorphic part of a LLL wave function. 

Specially, the anisotropic counterpart of the 
rotationally-invariant non-interacting basis states 

M^(7 = 0), s[^(7 = 0), (77) 

are 

M^(7), Ml)- (78) 

Therefore, if we could expand the many-body state 
\E'(7 = 0,ri, . . . ,rjv) in the non-interacting basis states 
M^{Q) (HI ors[^(0) (m, 

^cWX^(O), ^c^sI^(O), (79) 



its anisotropic counterpart ^(7, ri, . . . , tat) would be 

Y.'^fM.ii), Y.'^Ml)- (80) 

Here we find the coefficients and c^^^ are independent 
of 7. 

The 7-independence of the coefficients results in 7- 
independence of the occupation number for each orbital 

This also suggests that the entanglement spectrum^ 
which encodes the topological properties of the state, is 
also invariant, if the appropriate cut in the space of the 
total angular momentum associated with guiding cen- 
ter is chosen. We see that topological properties of a 
FQH state are built into the coefficients (relative weight), 
which are manifested, e.g., in the entanglement spectrum 
calculation. On the other hand, geometrical properties 
of the anisotropic FQH states are encoded into the de- 
formation of the non-interacting basis. Nevertheless, ge- 
ometrical information can be revealed by exploring the 
properties of the family of variational wave functions and 
the corresponding variational energies; this point will be 
illustrated by examples in Sec. HVl 

C. Anisotropic FQH states 

In this subsection, let us consider some prominent vari- 
ational functions such as the Laughlin wave function, 
Moore- Read state and so on. 

The V = l/q Laughlin wave function is given by 

*I(7 = 0) = n(^« - z,yJl^m{Q.r,). (81) 

i<j i 

and its anisotropic counterpart is given by 

Kii) = n^™(^'^^) n - + ^0(5^. - 5^.)]' 1- 

i i<j 

(82) 

The anisotropic Laughlin states ^^(7) are the (minimum 
total £(7)) zero-energy ground state of the Haldane's 
pseudopotential Hamiltonian i^int({Kn}, 7) with 

Vm<q > and Vm>q = 0. (83) 

To prove this, we may consider the diagonalization of 
the projected interaction Hamiltonian matrix and note 
the 7-independence of the matrix elements 

{BlMHint{{Vm},l)\Mh)), or 

(>^^(7)|i^int({K^},7)l^1p'(7)), 

which results from the 7-independence of 
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For two particles, the anisotropic counterpart of the 
V = \/q Laughhn state 



(84) 



is given by 



= V^oo(7,n)V^oo(7,^2) [^1-^2 + ^0(5^1 -5^2)]^!, 

= V^oo(7,^i)V^oo(7,^2)(^i2 + 22;^S^,2)n, 

= V^oo(7,^i)V^oo(7,^2)W(?(^i2,22;o), (85) 

with Z12 = z\—Z2 and 9^12 = \ (^^i — d^^. Among these 
states, the straightforward (anisotropic) generahzation of 
the = 1/2 bosonic Laughhn state is ^02(7, ''"1, ''"2), in 
which the double zero of the undeformed state at z\ = Z2 
is spht into zeroes at zi — Z2 = izV^zo- While for the 
u = 1/3 fermionic Laughhn state ^03(7, ''"1, ''"2), the triple 
zero at zi = Z2 is split into a single zero at zi = 2^2, and 
displaced zeroes at zi — Z2 = ±z\/6zo. This suggests the 
splitting in the pattern of zeros can be used as a tool 
to study the anisotropic deformation of a quantum Hall 
wave function. 

The other FQH state, such as the Moore-Read stated, 
the Read-Rezayi state^, the Haldane-Rezayi state^^, etc., 
could be generalized in this way, as long as the corre- 
sponding wave functions could be expanded in the non- 
interacting LLL basis. In addition, we also conclude that 
the anisotropic Zk parafermion states^ are also the (min- 
imum total generalized angular momentum) zero-energy 
ground state of some special {k + l)-body interaction. 
This kind of interaction could be expressed in terms of 
the generalized (/c-h l)-body pseudopotentials^ since the 
Zfe parafermion states are the (minimum total angular 
momentum) zero-energy ground state of {k + l)-body 
short-range interaction^'^^. 

For example, the v = 1/q anisotropic Moore- Read 
states are given by 



Y[ [z, - Zj + zlid,, - d,^)Y 1, (86) 

i<j 

where Pf means the Pfaffian, the square root of the de- 
terminant. Among these states, the anisotropic u = 1 
bosonic Moore- Read states ^^^(7) (minimum to- 

tal £(7)) zero-energy ground state of the anisotropic 3- 
body pseudopotential Hamiltonian, 



i^int(7) = Voho{-f), Vo > 0, 
where the projection operator ho{j) is given by 

00 

^ ^ ^(mi,m2,m3|M)^ 

M=0 mim2rn3m4m5me 

, 7725, '^g) ^ 9 rni 9^2^^3 91716 91715 977, 



(87) 




FIG. 3: (Color online) The density profiles (in units of 
l/(27r^^)) of the = 1/3 anisotropic Laughlin state with 
7 = 1/2 (a) and 7 = (b) for an AT = 10 system. Note 
the yellow bulk whose value approximates 1/3. 



with the 3-body basis state |mi , m2, ^3)^ = | V^o,mi (7)) ^ 
IV^o,m2(7)) ^ 1^^0,7713(7)) and 



|M), 



• <2 + <3 



M 



|0,0,0). 



And the anisotropic u = 1/2 fermionic Moore- Read 
states ^^^(7) is the (minimum £(7)) zero-energy 
ground state of the following anisotropic 3-body pseu- 
dopotential Hamiltonian, 

where the projection operator hi{j) is given by 



E E 



^(mi,m2,m3|3,M - 3)^ 



M = 771177127713777477757776 



^(3, M - 3|m4, m3,m6)^gl^^gl^^gl^^gme9m59m4- 
with the three-body state |3, M — 3)^ being 

|3,M-3)^ = BM{a\^-a\^){a\^-a\^){a\^-a\^) 
+4^ + 4 J |0,0,0)^ 
and Bm being the appropriate normalization factor. 



IV. NUMERICAL RESULTS FOR THE 
ANISOTROPIC FQH STATES 

In this section we numerically study some properties of 
the the anisotropic Laughlin states ^^(7) and its applica- 
bility in a system with dipole-dipole interaction. We first 
demonstrate the deformation of the FQH droplet in den- 
sity profile. The anisotropy leads also to the deformation 
of the correlation hole, which can be understood by the 
split of the third-order zero to three adjacent first-order 
zeros. As a trivial example we show that the isotropic 
Laughlin states, among the family of generalized states, 
is the ground state in the variational sense for isotropic 
hard-core interaction. On the other hand, the anisotropic 
dipole-dipole interaction picks a variational ground state 
with a finite 7 parameter as expected. 



8 



10 



1% 



10 -10 




FIG. 4: (Color online) The pair correlation function for u — 
1/3 Laughlin state with 7 = 1/2 (a) and 7 = (b) of an 
N — 10 system. 



A. Density profiles 

The density profile for the anisotropic Laughlin state 
is most easily calculated using the Jack polynomials^i^^. 
Explicitly, we can write 



(89) 



for a finite number of particles. In Fig. [3l we plot the 
density profile gij^z) with 7 = 1/2 for an A/" = 10 
fermionic system at u = 1/3 and compare with that of 
the isotropic Laughlin wave function. Roughly speaking, 
in the anisotropic case the FQH droplet is stretched along 
the ^-direction, while it maintains its value around 1/3, 
in units of l/(27r£^) in the bulk (indicated by the yellow 
color in the color plot). 



B. Pair correlation function 

The density-density correlation function represents the 
conditional probability of find one particle at r when an- 
other is simultaneously at r^ For the anisotropic Laugh- 
lin state ^1(7), the density-density correlation function 
is defined by 

G(2)(7,r,r') = {r)¥ {r')i>{r')i>{r)\^l{j). 

Without loss of generality, we consider the pair correla- 
tion function with fixed at the origin. 



^(7,r)^(7,0)' 



(90) 



In Fig. m we compare g{'y,r) of slu N = 10 Laugh- 
lin state at u = 1/3 with 7 = 1/2 and 0. The two 
cases are clearly distinguishable, as the correlation hole 
for 7 = 1/2 is stretched non-monotonically along the y- 
direction, along which the density profile is stretched. 
The center of the hole is strictly zero but two other min- 
ima developed in the ^-direction are not. To explore 
more details we analyze ^(7,2:, 7/) as a function of x (or 
y) along the y-ax.is (x-axis) in Fig. [5l The comparison 
shows that the asymptotic behavior of ^(7 7^ 0, r -> 0) is 



1.6 



1.2 



0.8 



0.4 



^(7,^,0) 



.X 10" 



(a) 



-15 -10 




-5 5 

x{i) or y{i) 



FIG. 5: (Color online) The pair correlation function ^(7, 0, y), 
g{j, X, 0) with 7=1/2 and ^(0, x, 0) of an iV = 10 anisotropic 
Laughlin state. The inset (a) shows the asymptotic (green 
solid) line of g{0,x,0), which is proportional to . And the 
inset (b) shows the asymptotic (green solid) line of ^(7,x,0) 
(^(7,0,2/)), which is proportional to (y^)- 



very different from that of ^(7 = 0,r ^ 0). The differ- 
ence roots in the toy model for two particles we studied in 
Subsec. UlI C[ from which we expect that g{j 7^ 0, r ^ 0) 
vanishes as |rp, while ^(7 = 0,r ^ 0) as |r|^, as indi- 
cated by Eqs. (|84|) and (|H5|). The insets (a) and (b) in 
Fig. [5] confirm the asymptotic behavior. 

We emphasize that in the anisotropic case the pair cor- 
relation function is isotropic for small enough |r|, reflect- 
ing the first-order zero in the fermionic wave function en- 
forced by the Pauli exclusion principle. The two-particle 
wave function can be regarded as the asymptotic wave 
function when two particles are significantly closer than 
their distance to any other particles. For \r\ comparable 
to I7I, the pair correlation function becomes anisotropic. 
The observed anisotropy at short distances encodes the 
geometrical deformation of the Laughlin state. When two 
particles are close to each other, each particle sees q zeros 
(for the Laughlin state at z/ = 1/q)] the spatial spread of 
them reflects the extent of the deformation. Topological 
properties, manifested in the isotropic case by the gth- 
order zero, can be identified when one look not too close 
i\r\ » l7l). 



C. Variational energy 

Given the family of Laughlin states ^^(7), we test 
the variational principle on a trivial Hamiltonian with 
the isotropic hard-core interaction, which renders the 
isotropic Laughlin state as its exact zero-energy ground 
state, and a Hamiltonian with anisotropic dipolar-dipolar 
interaction. In both cases the expectation value of the 
Hamiltonian develops a minimum, which may be iden- 
tified as the variational ground state energy. The mini- 
mum occurs at 7 = for the isotropic interaction and a 
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FIG. 6: (Color online) The hard-core interaction energy per 



particle 4^(7) with respect to the distort parameter 7 for pjQ 7. ((^^^^^ oxiXy^q) The dipole-dipole interaction energy 



iV = 3 
parison 



system. Here we also plot 7 (solid line) for com- 



nonzero 7 for the anisotropic interaction. 

Because the Laughhn wave function is the zero-energy 
ground state of the (isotropic) hard-core interaction, 
7 = is naturally the minimum for the energy expec- 
tation value. The variational energy per particle is ex- 
pected to increase as /3|7p for small I7I, as can be quickly 
understood, e.g., from the expansion of ^^(7) with re- 
spect to small I7I, 



1^* 
2 ' 



Y..dl + ..\ n(o).(9i) 



Here /3, like stiffness, quantifies the energy cost of the 
fluctuation of metric. In Fig. [6] we plot the hard-core in- 
teraction energy per particle for the anisotropic fermionic 
Laughlin state ^^"^(7) 



(vl/i(7)|gint({V„^},7 = 0)l^i(7)) 



,(92) 



with Vm = ViSm,i' For the set of finite-size systems with 
N = 3-8, we confirm that the minimum of the hard-core 
interaction energy occurs identically at 7 = 0. There- 
fore, the (isotropic) Laughlin wave function is indeed the 
optimal state for the hard-core interaction. 

We now turn to the dipole-dipole interaction, which 
some of us studied in Ref. 9 by exact diagonalization. For 
the dipole-dipole interaction with the dipole moments be- 
ing polarized in the x-z plane, the potential in the x-y 
plane is given by 



V{0,x,y) 



Cd 



(X2+ ^2+^2)5/2 

{x^ + + - 3(^ cos 6> + X sin of} , 



where Cd is the interaction strength and the motion of all 
the particles along the z-axis is frozen to the ground state 
of the axial harmonic oscillator 7r~^^^d~^^'^e~^ with 
d being the axial oscillator length. Here the tilt angle d is 
introduced to tune the dipole-dipole interaction such that 
V{B^ z) is isotropic (anisotropic) on x-y plane for 6> = 



per particle ^^^\0 = 30°, 7) with respect to the distort pa- 



rameter 7 for N 
d is set as 0.01-£. 



system. Here axial oscillator length 



{6^0). The dipole-dipole interaction energy per particle 
in the anisotropic fermionic Laughlin state ^^"^(7) is 



N 



(93) 



In Fig. El we plot S^J^\o = 30°, 7) with respect to 7 
for finite-size system with TV = 3 — 8. We found that 
the lowest ground state energy of the dipole-dipole inter- 
action occurs for the anisotropic Laughlin state with a 
nonzero 7. The further quantitative comparison between 
the model anisotropic quantum Hall state and the exact 
ground state of the anisotropic dipole-dipole interaction 
will be given elsewhere. 



V. SUMMARY 



In this work we have explicitly constructed families of 
anisotropic fractional quantum Hall states, which are ex- 
act zero energy ground states of appropriate anisotropic 
short-range two- or multi-particle interactions. These 
families generalize the celebrated Laughlin, Moore-Read 
and Read-Rezayi states. Each family is parameterized 
by a single geometric parameter that describes the dis- 
tortion of correlation hole in the density-density corre- 
lation functions. These states thus explicitly illustrate 
the existence of geometric degree of freedom in fractional 
quantum Hall effect, as recently demonstrated in Ref. 0. 
Application of these states in studies of systems with 
realistic anisotropic interactions, like dipole-dipole inter- 
action, will be pursued in the near future and reported 
elsewhere. 
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